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ABSTRACT 

We analyze the C4 catalog of galaxy clusters from the Sloan Digital Sky 
Survey (SDSS) to investigate the axis-ratio distribution of the projected two 
dimensional cluster profiles. We consider only those objects in the catalog whose 
virial mass is close to IO^'^H'^Mq, with member galaxies within the scale radius 
1000 kpc. The total number of such objects turns out to be 336. We also 
derive a theoretical distribution by incorporating the effect of projection onto 
the sky into the analytic formalism proposed recently by Lee, Jing, & Suto. 
The theoretical distribution of the cluster axis-ratios is shown to depend on the 
amplitude of the linear power spectrum (ag) as well as the density parameter 
(fim). Finally, fitting the observational data to the analytic distribution with 
fljn and CTg as two adjustable free parameters, we find the best-fitting value of 
as = (LOl ± 0.09) (1^^/0.6) (°-°^±°-°2)+o-if2™. it is a new (Xs-^m relation, different 
from the previous one derived from the local abundance of X-ray clusters. We 
expect that the axis-ratio distribution of galaxy clusters, if combined with the 
local abundance of clusters, may put simultaneous constraints on erg and flrn- 

Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

The standard theory of structure formation based on the cold dark matter (CDM) 
paradigm explains that the large scale structures of the universe like galaxies, groups and 
clusters of galaxies are originated from the primordial fluctuations of the matter density field 
through gravitational instability. In this scenario, the observed properties of the large scale 
structures can be used to extract crucial information on the initial conditions of the early 
universe. 
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Among the various properties of the large scale structures, it is the abundance of galaxy 
clusters that has attracted most cosmological attentions so far. Since the clusters are rela- 
tively young still in the quasi-linear regime, their abundance can be inferred from the linear 
theory (e.g., Press & Schechter 1974). Furthermore, being rare, the evolution of galaxy clus- 
ters depends sensitively on the background cosmology, especially on the amplitude of the 
matter power spectrum on scale of 8h^^ Mpc (cts) and the density parameter (flm)- It has 
been shown that by comparing the linear theory prediction with the cluster abundance from 
X-ray or Sunyaev-Zel'dovich (SZ) effect observations one can constrain erg and (White 
et al. 1993; Bond & Myers 1996; Eke et al. 1996; Pen 1998; Henry 2000; Fan & Chiueh 
2001; Pierpaoh etal. 2001; Seljak 2002). For instance. Eke et al. (1996) found a relation of 
as — (0.52 ± 0.04)Q^'^'^~^^'^^^"^ (for a flat universe with non-zero cosmological constant. A) 
by comparing the Press-Schechter prediction with the local abundance of the X-ray clusters 
complied by Henry & Arnaud (1991). 

To break the degeneracy between as and 0,^, however, it is necessary to find another 
relation obtained from another observable property of galaxy clusters. One candidate is the 
shape of galaxy clusters. It has long been known observationally as well as numerically that 
in general the shapes of galaxy clusters are not so much spherical as elliptical (Schipper 
& King 1978; Binggelh 1982; West et al. 1989; Jing & Suto 2002). Furthermore, it was 
shown by recent numerical simulations that the ellipticity distribution of galaxy clusters 
depend strongly on the background cosmology (Jing & Suto 2002; Hopkins etal. 2005). This 
numerical finding implies that it may be possible to find a new independent relation between 

and as with the ellipticity distribution of galaxy clusters. 

Most of the previous approaches to the cluster ellipticity distribution were numerical. In 
order to use the cluster ellipticity distribution as a cosmological probe, however, it is highly 
desirable to have an analytic model derived from physical principles. Although Bardeen et 
al. (1986) made a first analytic attempt in the frame of the density peak formalism, their 
model was applicable only in the asymptotic limit of low ellipticity. 

Very recently, Lee et al. (2005, hereafter, LJS05) derived a new analytic expression for 
the halo axis-ratio distribution by combining the density peak formalism and the Zel'dovich 
approximation. They demonstrated that the analytic model produces the characteristic 
behaviors of the halo axis-ratio distribution found in simulations, which gives a hope that 
it may provide a theoretical footing for the use of the cluster ellipticity distribution as a 
cosmological probe. Our goal here is to find a new functional form of the aQ—Q^ relation, 
independent of the previous one derived from the cluster abundance, by comparing the 
LJS05 analytic prediction to the axis-ratio distribution of galaxy clusters determined from 
the recent observational catalog. 
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The plan of this paper is as follows: In §2, we provide a brief description of the ob- 
servational data. In §3, wc explain how to derive an analytic axis-ratio distribution of two 
dimensional projected profiles of galaxy clusters from the LJS05 formalism. In §4, we find 
a new functional form of relation (Tg—Qm relation through fitting the model to the observa- 
tional data. In §5, we discuss the success and the limitation of our work, and draw a final 
conclusion. 



2. DATA FROM SDSS-C4 CATALOG 

The C4 catalog contains 748 objects (Miller et al. 2005) identified in the Second Data 
Release (DR2) of the Sloan Digital Sky Survey (SDSS, Strauss et al. 2002). The sky coverage 
and the rcdshift range of the catalog arc ~ 2600 dcg^ and 0.02 < z < 0.17, respectively. The 
majority of the objects in this catalog are clusters of galaxies with the mean membership of 
36 galaxies. In the catalog, however, are also included smaller groups with the membership 
of at most 10 galaxies as well as larger structures with the membership of over 100 clusters. 
The mass range of the catalog is thus quite broad from 10^^/i~^Moto lO^^/i^^M© {h: the 
dimensionless Hubble parameter). 

The SDSS-C4 catalog provides various spectroscopic properties of galaxy clusters, among 
which the following informations are used for our study: the axis-ratio of a projected two- 
dimensional image (q), the virial mass (M), and the redshift (z). As below, we briefly 
describe how these quantities were measured. For the detailed descriptions of the C4 catalog 
and the measurement of M and q, refer to the web-page: 

http://www.ctio.noao.edu/~chrism/current/research/C4/dr2/ SDSS_PARAMETERS.html. 

The axis ratio, q, is derived from an elliptical fit to a projected two dimensional image 
of a cluster. The elliptical fit was conducted by the diagonalization of the covariance of the 
positions of all galaxies within 1000 kpc with respect to the cluster center. The center of a 
cluster is also measured using all galaxies within 1000 kpc. 

The virial mass, M, of each cluster was computed (assuming the isothermality and the 
spherical shape) as 

M = ^afr, (1) 

where ai is the line-of-sight velocity dispersion and is the deprojected three dimensional 
virial radius of a cluster (Carlberg et al. 1996). The value of r-^ is related to the ringwise 
projected harmonic mean radius, Rh (see eq. 3 in, Carlberg et al. 1996), as = T:Rh/2 
(Limber & Mathews 1960). The value of Rh depends on which scale radius is used in selecting 
the member galaxies. In the catalog, there are five different criteria used to measure the 
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value of M: w = 500, 1000, 1500, 2000, and 2500 kpc. 

In this paper, consistently, we use the criterion w — 1000 kpc. Using the fixed scale 
radius w — 1000 kpc does not mean that the virial radius of a cluster is fixed to be 1000 kpc 
for all clusters. What it means is that the ringwise projected harmonic mean radius, Rh, 
was determined using those galaxies within the scale radius w — 1000 kpc from the cluster 
center. 

But, note that using the fixed scale radius oi w — 1000 kpc may systematically affect 
the determination of the cluster mass, M. Since the catalog contains various objects in the 
wide mass range, using the fixed scale radius oi w — 1000 can cause either underestimate 
or overestimate of the vahic of M, depending on how massive an object is. To reduce 
this kind of possible systematic effect, we focus only on the narrow cluster mass range of 
< M < We find total 336 such clusters in the catalog with the 

median redshift of Zm. = 0.08. 



3. PHYSICAL ANALYSIS 

3.1. Overview of the Analytic Model 

The LJS05 analytic model for the axis-ratio distribution of dark halos adopts two clas- 
sical theories as its background: the Zel'dovich approximation (Zeldovich 1970) and the 
density peak formalism (Bardeen et al. 1986). 

By applying the Zel'dovich approximation, LJS05 assumed that the three principal axes 
of a triaxial halo, {a, b, c} (with a < b < c) are related to the three eigenvalues {Ai, A2, A3} 
(with Al > A2 > A3) of the linear deformation tensor, dij = didj(j) where (p is the perturbation 
potential: 

a oc y/1 - L>+Ai, b oc y/1 - D+X2, C oc y/1 - D+X3, (2) 

where — D+{z) is the growth rate of the hnear density field. Here, the sum of the three 
eigenvalues equals the linear density contrast, 5 = Ap/p (p: the mean density): A1-I-A2+A3 = 
S. 

By employing the density peak formalism, LJS05 also assumed that a region in the linear 
density field will condense out a virialized halo if it satisfies the following two conditions: 

6 = 6, = 6,oD40)/D4z), (3) 
A3 > Ao = X,oD+{0)/D+{z), (4) 

where 5c and Ac represent the threshold values of 5 and A, respectively, with = Sc{z — 0) 
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and Aco = Xc{z — 0). The value of Sco can be theoretically evaluated with the help of the 
top-hat spherical collapse model. For instance, 5co ~ 1.686 for a flat universe with closure 
density (e.g., Kitayama & Suto 1996). 

As for the value of Aco, it is precisely in the original density peak formalism. Instead 
of using this theoretical value, however, LJS05 used a positive value of Ac = 0.37 which 
was determined empirically through fitting. Their claim was that the break-down of the 
linear theory in the highly nonlinear regime is responsible for the deviation of Aco from the 
theoretical value of in practice. 

Here, we use the theoretical value of Aco = rather than the empirical value of Aco = 0.37 
for the following reason. The value of Aco = 0.37 was found empirically by comparing the 
analytic model with the numerical result for the concordance cosmology only. Remember 
that our purpose is to use the analytic distribution as a probe of cosmology without having 
any bias. Hence, we set Aco at the original theoretical value in the rest of this paper. 



Defining two real variables /ii and 112 as 

b 

1^1 = 



a 
c 



(5) 



LJS05 derived analytically the joint probability density distribution of /ii and 112 which 
depends on the cluster mass M and the formation epoch Zf: 



, 3375^2 
A ^ ^ -exp 



M 



5Sl 



2< 



+ 



155,(Ai + A2) 15{XI + A1A2 + Ai) 



2< 



X (2Ai + A2 - 6c) (Ai - A2) (Ai + 2A2 - 5c) 

{dXid\2) 



X 



e(i57^)ei4-(A. + A.) 



2< 



(6) 



where Df = D{zf), a is the standard variation of the linear density field, and A is the 
normalization factor which satisfies 



^ j p{lJ'i,IJ'2;M,Zf)diJ,idfj.2 = 1. 



(7) 



The relation between {/ii,/i2} and {Ai, A2} arc given as 

l + {Df6c-2)fil + fj,i 



Ai = 
A2 = 



l + iDfSc-2)ixl + iJ,l 



Dfii^l + lil + l) 



(8) 
(9) 
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and the Jacobian |(9Ai9A2)/(9/xi9//2)| was found to be 



{d\id\2 

((9/ii5/i2 



1)= 



(10) 



Equation (6) gives the probabihty distribution of the two axis-ratios of a triaxial galaxy 
cluster with mass M formed at 2;/. In practice, however, what we measure is not the 
formation epoch, Zf but the observation epoch, z, i.e., the redshift at which a cluster is 
observed. Hence, LJS05 rederived an analytic expression for p{iii, 112', M; z) from equation 
(6) as 

f°° dvfizf- 2M z) 
p{l^i,l^2;M;z) = J dzf ^-^^ ' ' p{i^i,i^2]2M;zf), (11) 

where dpf/dzf represents the formation epoch distribution of galaxy clusters. LJS05 used 
the fitting formula given by Kitayama & Suto (1996) which is an approximation to the 
analytic expression found by Lacey & Cole (1994). 



3.2. Projection onto the Sky 

Although equations (6)-(ll) provides an analytic expression for the axis-ratio distribu- 
tion of the three dimensional shape of a triaxial cluster with mass M observed at z, it is 
hard to compare it directly with the observational data. In observation, as we noted in §2, 
what is measured is only a two dimensional projected profile of a galaxy cluster. Therefore, 
for a proper comparison with observational data, it is necessary to incorporate the effect of 
projection along the line of sight into the analytic model. 

To incorporate into equation (11) the projection effect onto the plane of the sky, we 
follow Binney (1985): Let {9, (p) be the usual polar coordinates of the hne-of-sight vector in 
the principal axes of a triaxial cluster, and let q be the axis-ratio of a two dimensional cluster 
profile projected along the line of sight direction. Then, the probabihty density distribution 
of q can be obtained by performing an integration of equation (6) over //i, //2, ^, as 

p{q;M,z)^— d(f) sinOde d/ii dii2 Soiq - q'{0,(f), /ii, 112)] p{lJ'i, IJ'2; M, z), {12) 
47rJo Jo Jn2 Jo 



where 



(13) 
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and 

008^9 f . 2^ cos^0\ sin^^ ,^ ^, 

U = —^[sm'9 + —^\+—^, (14) 

W = cos^sin20 ( 1 - ^ J ^, (15) 

/ sin^ 6 r, ,\ 1 
V . (^_^ + oos^,j^, (16) 

Through equations (6)- (16), one can evaluate analytically the probabihty density that the 
projected profile of a galaxy cluster of mass M is observed to have an axis-ratio of q at 
redshift z. This theoretical axis-ratio distribution varies with the background cosmology 
since it depends on the initial power spectrum, P{k). In this paper, we adopt the following 
approximation formula for P{k) (Bardeen et al. 1986): 

>n 2 



P{k) oc k 



ln(l + 2Mq) 
2.34? 



[1 + 3.89g + (16.1g)' + (5.46)-^ + (6.71g) V'^', (17) 



where q = k/Th ^ Mpc. Here, the shape factor, F, is related to Jl^, h, and the baryon 
density parameter, Qi, (Peacock & Dodds 1994; Sugiyama 1995) as 

where Tq is the temperature of the cosmic microwave background radiation. 

Fig. 2 plots equation (12) on cluster scale, Mg, (see §4 for a definition of Mg) at present 
redshift {z = 0) for four different cosmological models: ACDM (solid line); OCDM (dotted); 
SCDM (dashed); rCDM (White et al. 1995, long-dashed). The cosmological parameters 
used to characterize each model are listed in Table 1. 

As one can see, equation (12) depends sensitively on the background cosmology. More 
elliptical clusters are predicted by the ACDM and OCDM models with high ag and low Qm 
than by the SCDM and rCDM models with low as and high flm- It can be understood by 
the following logic: the SCDM and rCDM models predict more small scale powers, and thus 
should result in rounder shapes of galaxy clusters in the end. 

We compare equation (12) with the observational distribution from the 336 SDSS-C4 
clusters, assuming the ACDM. Figure 2 plots together the analytic (solid line) and the 
observational (dots) axis-ratio results for comparison. The errors of the observational points 
are Poissonian. The values of M and z in the analytic axis-ratio distribution are set to be 
Ms (see §4.1) and Zmed — 0.08, to be consistent with the observational data. As one can see, 
the overall agreement between the two is quite good. 
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4. APPLICATION OF THEORY TO OBSERVATIONS 

4.1. Characteristic Cluster Mass 

The analytic axis-ratio distribution (eq.[12]) depends on the cluster mass M and redshift 
z. Therefore, to compare it with the observational data, one has to first specify the values 
of M and z. As for z, to be consistent with the observational data, we set it at the median 
redshift of the selected 336 clusters, z^ — 0.08, given that the redshift range of the selected 
clusters is quite narrow. 

As for M, there is a subtle point that has to be taken into consideration. The virial mass, 
M, of each cluster in the C4 catalog was computed under the assumption of a concordance 
cosmology with Q,n = 0.3, <ts = 0.9, h — 0.7. In other words, the cluster mass M given 
in the catalog is biased toward the concordance cosmology, which has to be avoided in the 
parameter estimation. Thus, we set M at the characteristic mass scale, Mg, defined as 

Ms = —pRl (19) 

where Rs = 8/i~^Mpc is the top-hat spherical radius on which scale the rms density fiuc- 
tuation is observed to be very close to unity, and p = 2.78 • 10^^ flmh'^ MqMpc~^ . Since the 
galaxy clusters are the largest collapsed objects in the universe, the characteristic mass scale 
Mg is believed to represent the typical mass of galaxy clusters. Note that this characteristic 
cluster mass, Mg, by its definition (eq.[19]), depends on the value of flm- Therefore, putting 
M — Mg into equation (12), will optimize the parameter estimation without bearing any 
bias toward a certain cosmology. 



4.2. Constraint on and (7g for a Flat Universe 

The analytic distribution (12) depends not only on Qm and (jg but also on Q\, h, and F 
as well. Here, we focus mainly on constraining fl^n and (jg, assuming that the other param- 
eters are already known as priors. We use the following priors from the CMB observation 
(e.g., Lange et al. 2001), the HST Key Project (Preedman et al. 2001) and the big bang 
nucleosynthesis (Olive etal. 2000): Qa = 1 — (a flat universe); h — 0.7; Qb — 0.044. Once 
the values of h and are given, the shape factor, F, will depend only on the value of 
(see eq.[18]). 

Now that all parameters except for fi^ and erg in equation (12) are prescribed, we can 
fit the observational data points obtained in §2 to equation (12) with adjusting and (7g 
in ranges of (0, 1] and (0, 2], respectively. We pick as the best-fit values of and erg those 
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that minimize : 

" h ' ^ ^ 

where A^p is the total number of observational data points, {qi,Pob{(li)}, and aob{qi) is the 
observational standard variation of qi. 

The single best-fit values, however, are unreliable given the uncertainty in due to 
the strong correlation between Vim and ag. A more general best-fit is defined as the full 
range of the outer limit of one standard deviation contour. Figure 3 plots the one standard 
deviation contour as solid line. We find that this contour is well approximated by the 
relation, ag = (1.01 ± 0.09)((]„/0.6)(°-°^±°-°2)+°-^"-. This approximation formula is also 
plotted as dot-dashed line in Figure 3. The previous relation of Us = (0.52±0.04)Q~°-^^+°-^^^'" 
obtained by Eke et al. (1996) from the cluster abundance is also plotted for comparison as 
dashed line. The shaded area surrounded by the solid and the dashed lines in Figure 3 
represents the simultaneously constrained best-fit values of and cxg: = 0.31 ± 0.07 

and = 0.94 ± 0.07, which is in good agreement with the results from the first year 

WMAP (Wilkinson Microwave Anisotropy Probe) measurements (Spergel et al. 2003). 



5. DISCUSSION AND CONCLUSIONS 



We have found a new as—^m relation by comparing the LJS05 analytic model for the 
axis-ratio distribution of galaxy clusters with the observational data from the SDSS-C4 
catalog. 

The success of our final result, however, is subject to a couple of caveats. The first 
caveat comes from our assumption that the principal axes of the SDSS-C4 clusters derived 
from the elliptical fits to the galaxy distributions are directly comparable to that of the 
analytically derived shapes of dark halos. This is an obvious simphfication of the reality. 

Notwithstanding, the following observational and numerical clues should provide a jus- 
tification for this simplified assumption. Many observations revealed that the spatial distri- 
bution of cluster galaxies and the major axes of their host clusters are arranged in a coUinear 
way (West & Blakeslee 2000; Phonis & Basilakos 2002; Plionis et al. 2003; Pereira & Kuhn 
2004). Recent numerical simulations also demonstrated that the dark halo substructures are 
preferentially located along the major axes of their host halos (Kang et al. 2005; Libeskind 
et al. 2005; Zentner et al. 2005; Lee & Kang 2006). Very recently, Lee & Kang (2006) showed 
by N-body simulations that the triaxial shapes of dark matter halos can be reconstructed 
from the the anisotropic spatial distribution of their substructures. Given these empirical 
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findings, one may expect that the cluster principal axis derived from the elliptical fit to the 
galaxy distribution is comparable to the analytical one derived from the the smooth dark 
matter distribution. 

The second caveat lies in the hmited vahdity of the analytic model. As shown by LJS05, 
the prediction of our analytic model on the cluster ellipticity-mass relation disagrees with the 
numerical finding. Our model predicts that the cluster ellipticity decreases with mass, while 
many different numerical simulations (Bullock 2002; Jing & Suto 2002; Kasun & Evrard 
2005; Hopkins etal. 2005; Bailin & Steinmetz 2005) found that the cluster ellipticity actually 
increases with mass. 

In spite of this limitation of the LJS05 analytic model, it is still true that the LJS05 
tested their model against N-body simulations and showed that it works quite well, repro- 
ducing the characteristic behaviors of the numerical result, on a single mass scale (or in 
a narrow mass range). In other words, although the LJS05 model fails in predicting the 
"change" of the axis-ratio distribution in the wide range of halo mass, it can still be used to 
predict the shape of the axis-ratio distribution of a dark halo with a given mass. 

As we consider here only those objects in the C4 catalog whose mass lies in a narrow 
range ~ 10^^h~^MQ for the comparison with the analytic model, the LJS05 model is expected 
to be a useful approximation. 

Definitely, however, it will be desirable and necessary to refine our model more realisti- 
cally before applying it to a wide range of mass of cosmic structures. Finally, we conclude 
that it is possible in principle to constrain both Q^n and as concurrently from the observa- 
tion of galaxy clusters alone, as the cluster shape and abundance distributions provide two 
different erg— il^ relations. The amplitude of the linear power spectrum obtained from the 
shape and abundance distribution of the galaxy clusters is consistent with the high value 
estimated from the recent WMAP measurements. 

We thank an anonymous referee for careful reading and many helpful suggestions that 
helped us greatly improve the original manuscript. We acknowledge stimulating discussions 
with R. Sheth, M. Im, C.Park, Y. Jing, Y. Choi, H. Hwang, J. J. Lee. This work is supported 
by the research grant No. ROl-2005-000-10610-0 from the Basic Research Program of the 
Korea Science and Engineering Foundation. 
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Table 1: The values of the cosmological parameters for each model 



model erg h T 

ACDM 0.3 0.7 OJ OJ 0.168 

OCDM 0.3 0.87 0.83 0.25 

SCDM 1 0.52 0.5 0.5 

rCDM 1 0.52 0.5 0.25 
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Fig. 1. — The axis-ratio distributions of the two dimensional projected profiles of galaxy 
clusters predicted by the analytic LJS05 model (eq.[12]) at present epoch for four different 
cosmological models: a ACDM model with Q„ = 0.3, Qa = 0.7,(T8 = 0.9,h = 0.7, F = 0.168 

(solid line); OCDM model with = 0.3, = 0, o-g = 0.87, h = 0.83, T = 0.25 (dotted); 
SCDM model with = 1, = 0, = 0.52, h = 0.5, T = 0.5 (dashed); rCDM model 
with = 1, Oa = 0, as = 0.52, h = 0.5, T = 0.25 (long-dashed). 
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Fig. 2. — The analytic axis-ratio distribution of the two dimensional projected profiles of 
galaxy clusters (solid line) and the observational data points arising from the SDSS-C4 
catalog (solid dots). The error bars are Poissonian. 
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Fig. 3. — The best-fit values of a^—Qm through x^-statistics. The solid and the dot-dashed 
lines represent the one standard deviation contour and its approximation. The dashed lines 
is the previous relation found from the local abundance of X-ray clusters. 



